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ABSTRACT 

We construct one dimensional steady-state models of accretion disks produced by the 
tidal disruption of a white dwarf (WD) by a neutron star (NS) or stellar mass black hole (BH). 
At radii r < 10 8 5 - 10 9 cm the midplane density and temperature are sufficiently high to burn 
the initial white dwarf material into increasingly heavier elements (e.g. Mg, Si, S, Ca, Fe, 
and Ni) at sequentially smaller radii. When the energy released by nuclear reactions is com- 
parable to that released gravitationally, we term the disk a nuclear-dominated accretion flow 
(NuDAF). At small radii < 10 7 cm iron photo-disintegrates into helium and then free nuclei, 
and in the very innermost disk cooling by neutrinos may be efficient. At the high accretion 
rates of relevance ~ 1CT 4 - 0.1 M & s _1 , most of the disk is radiatively inefficient and prone to 
outflows powered by viscous dissipation and nuclear burning. Outflow properties are calcu- 
lated by requiring that material in the midplane be marginally bound (Bernoulli constant < 0), 
due (in part) to cooling by matter escaping the disk. For reasonable assumptions regarding the 
properties of disk winds, we show that a significant fraction (> 50 - 80%) of the total WD 
mass is unbound. The composition of the ejecta is predominantly O, C, Si, Mg, Ne, Fe, and S 
[He, C, Si, S, Ar, and Fe], in the case of C-O [pure He] WDs, respectively, along with a small 
quantity ~ 10~ 3 - 1O~ 2 M0 of radioactive 56 Ni and, potentially, a trace amount of hydrogen. 
Depending on the pressure dependence of wind cooling, we find that the disk may be ther- 
mally unstable to nuclear burning, the likelihood of which increases for higher mass WDs. 
We use our results to evaluate possible electromagnetic counterparts of WD-NS/BH merg- 
ers, including optical transients powered by the radioactive decay of 56 Ni and radio transients 
powered by the interaction of the ejecta with the interstellar medium. We address whether 
recently discovered subluminous Type I supernovae result from WD-NS/BH mergers. Ulti- 
mately assessing the fate of these events requires global simulations of the disk evolution, 
which capture the complex interplay between nuclear burning, convection, and outflows. 
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1 INTRODUCTION 

Sensitive wide-field optical surveys are revolutionizing our under- 
standing of time-dependent astrophysical phenomena. New types 
of stellar explosions, that were once too faint or rare to be detected, 
are now routi nely discovered by efforts such as the Palomar Tran- 
sient Factory dLaw et alj|2009l ; iRau et alj[2009b and PanSTARRs 
dKaiseretalJl2002h . As our census of the transient universe ex- 
pands, it becomes increasingly important to evaluate the expected 
electromagnetic counterparts of known astrophysical events. 

Several Type I supernovae (SNe) have recently been discov- 
ered that are dimmer and/or more rapidly evolving than normal 



SNe la or Ib/c. Included among these is a class of peculiar la 
SNe characterized by lower luminosities an d lower ejecta velocities 
than n ormal la's ('SN 2002cx-like events' i llTet al.l2003l : |jha et al.1 
and which occur predominantly in star-forming galaxies 
dFolev et al]|2009h . SN 2008ha is an extreme example, with a peak 
brightness and rise time of only My ^ -14 and ~ 10 days, re- 
spectively dValenti et all 120091 : (Foley et al.l I2009L l2010h - Another 
recently identified class of subluminous transients are SNe lb with 
Ca-ri ch [but S, Si, and Fe-poor] spectra (e.g. SN 2005e; |Perets et al.l 
120101) . These events app ear to result from ejecta that have un- 
dergone helium burning dPerets et alj|2010l 1201 ll; IWaldman et all 
120101) and, unlike 2002cx-like events, are associated with an older 
stellar population located in the outskirts of their host galaxies. 
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Even more rapidl y-evolving SNe such as 2002bj dPoznanski et al.l 
l20ld) and 2010X dKasliwal et alj201Ch may be related to this class. 

Lacking hydrogen in their spectra, Type I SNe are generally 
thought to originate from compact progenitors such as white dwarfs 
(WDs) or massive stars without extended envelopes. Because com- 
pact stars lose most of their initial thermal energy to adiabatic 
expansion during the explosion, their emission must be powered 
by continued energy input from the radioactive decay of isotopes 
such as 56 Ni. The rapid evolution and low luminosities of events 
like 2002cx and 2005E thus require both a significantly lower 56 Ni 
yield, and total ejecta mass, than characterize normal SNe. 

It is economical to associate new classes of SNe with antici- 
pated variations of well-studied models, such as the core collapse 
of a massive star or the thermonuclear explosion of a WD. Events 
like 2002cx may, for instance, result from the pure deflagration of 
a Chadrasekhar or sub-Ch andrasekhar mass WD (e.g. branch et al.l 
l2004l;|Phillips et alj|2007l) . as opposed to th e 'delayed detonation ' 
models that best describe normal la SNe dNomotoetal.lll984h . 
They may alternatively result from we ak core collapse explosions 
dValenti et al.ll2009l : iMoriva et al.ll2010h . For Ca-rich lb SNe like 
2005E, a promising explanat ion is the detonation of a helium shell 
on the surface of a C/O WD dWooslev et al|l98j|llben & Tutukoyl 
19911; iLivne & Arnett|[l995l ; lBildsten et alj|2007l ; lshen et alj|20ld ; 
Wooslev & KaserJl20ld ). Th ough promising, none of these expla- 



nations is yet definitive (e.g. IWooslev & Kasenll20 1(X). The possi- 
bility thus remains that at least some of these events represent an 
entirely new type of stellar explosion. 

In this paper we examine the observable signatures of the 
merger of a WD with a binary neutron star (NS) or black hole (BH) 
companion. Four Galactic WD-NS binaries are currently known 
that will merge due to the emission of gravitational radiation within 
a Hubble time (see 5}2)- If mass transfer following Roche Lobe 
overflow is unstable, then the WD is tidally disrupted, producing 
a massive disk around the NS or BH. Past work has focused on 
the possibil ity that accretion powers a relativistic jet and a Gamma- 
Ray Burst dFrveret al.l 119991 iKing et alj|2007l) . Here we instead 
focus on the properties of the accretion flow on larger scales. We 
address two aspects of the problem neglected in previous studies: 

(1) the presence of [non-relativistic] outflows from the disk; and 

(2) the effects of nucleosynthesis on the thermodynamics and com- 
position of the disk and outflows. Although multidimensional nu- 
merical simulations are ultimately necessary to evaluate the evolu- 
tion and fate of these systems, the large range in spatial and tem- 
poral scales involved make such a study numerically challenging. 
Here, as a first step, we instead construct a simplified model of the 
disk and its outflows that we believe captures some of the essential 
physics. We use our results for the mass, velocity, and composition 
of the ejecta to quantify the associated electromagnetic counter- 
parts of WD-NS/BH mergers for a wide range of systems. 

This paper is organized as follows. In ^2]we discuss the con- 
ditions for unstable mass transfer in WD-NS and WD-BH binaries 
and motivate the initial conditions for our subsequent calculations. 
In 5j3] we describe one dimensional steady-state models of the ac- 
cretion disk and its outflows. In $4] we present our results for the 
disk structure 04.11 ) and outflow properties ( 94.2| l. In 94. 3 1 we as- 
sess the stability of our solutions and in 94.41 we discuss the ef- 
fects of convection. In ^5] we use our results to evaluate possible 
electromagnetic counterparts of WD-NS/WD-BH mergers, includ- 
ing subluminous Type I SNe ( 95. 2t and radio transients ( 95.3t . We 
conclude in §6\ 



2 UNSTABLE MASS TRANSFER AND DISK 
FORMATION 

A standard scenario for the formation of tight WD-NS/BH 
binaries invokes common envelope evolution of an initially 
wide binary, consisting of a neutron star or black hole and 
a int ermediate mass < 8 — IOMq ma in sequence companion 
(e.g. Ivan den Heuvel & Bonsdemlj 1 1984b . For a limited range of 
orbital periods, unstable Roche lobe overflow (RLOF) begins only 
after the secondary has left the main sequence. Depending on the 
evolutionary state of the stellar core when this occurs, the end re- 
sult is a NS or BH in orbit with a WD with either (from lowest 
to highest WD mass) a pure He, He-C-O ('hybrid'), C-O, or O-Ne 
composition. Close WD-NS/BH binaries can alternatively form di- 
rectly by collisions in dense stellar regions, such as th e centers of 
galaxies or globular clusters dSigurdsson & Reesll 19971) . 

After the WD-NS/BH binary is brought close together, 
the system continues to lose angular momentum on a longer 
timescale to gravitational wave emission. Angular momentum 
losses may be enhanced by the eccentricity indu ced by Kozai os - 
cillations if a tertiary companion is present (Tho mpson! l20ld). 
Of th e > 20 WD-NS binaries identified in our Galaxy dLorimerl 
2005), only four are sufficiently compact that they will merge 
in < 10'° yr (PSR J075 1+ 1807 - iKaspi et alj|200cj; J1757-5322 



- lEdwards& Bailed l200ll; Jl 14 1-6 545 - iLundgren etahl Il995t 
J1738-0333 - iFreire & WexbOld: see lO'Shaughnessv & Kiml2010l 
for a recent compilation). The WD masses in these systems are 
Mwd = 0.18, 0.67, 0.99, and 0.2M©, respect ively (e.g.lBailes et al.l 
2003). From the first three of these systems. iKim et al.ld2004l) esti- 
mate that the WD-NS merger rate in the Milky Way is 10~ 6 - 10~ 5 
yr , although correcting for pulsar beaming increases this rate 
by a factor of several. Population synthesis models predict some- 
what higher rates ~ 10~ 5 — 10~ 3 yr -1 , but with larger uncer- 



tainty ( Portegies Zwart & Yungelsonll 999: Tauri s & Sennelsl200d: 
iDavies et alj|2002h . No confirmed WD-BH systems are currently 
known. 

Once the orbital period decreases to < 1 minute, the WD un- 
dergoes RLOF onto the NS/BH companion. Approximating the 
WD as a T = 5/3 polytrope, its radius is 



10 4 



M 



WD 



0.1 M, 



-1/3 



1 - 



4/3 



V 



km, 



(1) 



where M C h = 1.4M e an d ji e is the mean molecular weight per elec- 
tron (jNauenberg 19721). The or bital separation at Roche contact is 
approximately (jEggleton 19831) 



0.6<r /3 +ln(l + q ,/3 ) 



0.49,? 



2/3 



(2) 



where q = M WD /M and M is the mass of the primary BH or NS. 

If orbital angular momentum is conserved, then mass trans- 
fer is unstable for q > 0.4 - 0.55 (e.g. lVerbunt & Rappaport] fl988t 
see lPaschalidis et alj|2009l . their Fig. 11). If strictly applicable, this 
criterion would limit tidal disruption to NS binaries with relatively 
massive > 0.5Mq C-O/O-Ne WDs, and might preclude disruption 
in WD-BH binaries altogether. Stable systems increase their orbital 
separation after mass transfer begins. The result is a long-lived ac- 
creting syst em, which may be observed as a n ultracompact X-ray 
binary (e.g. lVerbunt & van den Heuvelll 19951) . 

Conservative mass transfer is, however, unlikely. Tidal cou- 
pling during inspiral can transfer orbital angular momentum into 
spin. Furthermore, the accretion rate onto the NS/BH just af- 
ter RLOF is highly super-Eddington, in which case a com- 
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mon envelope may en gulf the system and outflows are likely 
(e.g. lOhsuga et alj2005l) . Depending on the amount of angular mo- 
mentum lost to winds or tides, unstable mass trans fer may occur in 
systems with even lower mass ratios q < 0.3 (e.g. lYungelson et all 
2002). In this case, tidal disruption may occur also for NS binaries 
with He or 'hybrid' He-C-O WD companion s, or even in the ca se of 
C-O/O-Ne WDs with low mass ~ 3Mq BHs jFrver et al.ll999h . Fu- 
ture numerical simulations are required to address t he challenging 
issue of stability in degene rate binary mergers (e.g. iGuerrero et al.l 
12004 iD'Souza et alfeOOfj) . 

If mass transfer is unstable, then the WD is tidally disrupted in 
just a few orbits. The WD material then circularizes and produces 
a disk around the central NS/BH with a total mass Mwd and char- 
acteristic radius R d ~ /?rlof(1 + <?)~ 3 - This disk accretes onto the 
NS/BH on the viscous timescale 



GM 



1/2 



140 s 



Get)" 



1.4M 



•1/2 



/ j? d \ 3l2 lH/Ri 
\10 9 cm/ \ 0.4 



(3) 



and at a characteristic rate 
M(R d ) ~ 



1O -M s-'f^)f^- 
u \0.1/\1.4M o 



1/2 



Mi 



0.7 Mr. 



UO 9 cm/ 



-3/2 



H/Rc 
0.4 



(4) 



where M is the NS/BH mass, a parametrizes the disk viscosity 
(Sj3}, and H is the scale-height of the disk, normalized to a char- 
acteristic (i)4}. Turbulence in the disk is driven by the magneto- 
rotational instability, which transports angular momentum at a rate 
corresponding to a ~ 0.01 - 0.1 (e.g. iDavis et a l. 2010). Because 
the mass of the disk is significant compared to that of the central 
NS or BH, the disk may also be gravitationally unstable, which 
likel y results in even more efficient angular momentum transport 
(e.g. lLaughlin & Bodenhei meill 19941) . In 94.4l we describe how the 
disk lifetime is altered if convection plays an important role in 
transporting angular momentum. 

Applying equations 10 and (|4j to the full range of unstable 
systems, we conclude that the majority of the WD mass accretes at 
a rate M(R A ) ~ 10~ 4 - 10~'Mq s -1 , depending primarily on M WD 
and a. More massive WDs accrete at a higher rate, because both 
M cc R d 312 oc R^ 2 F (l + q) 9 ' 2 and possibly a (if the disk is grav- 
itationally unstable) increase with q. At such high accretion rates 
(~ 10" - 10 14 times larger than the Eddington rate) the disk cannot 
cool through photon emission and is termed a radiatively-inefficient 
accretion flow (RIAF). 



3 ACCRETION DISK AND OUTFLOW MODEL 

In this section we present a one dimensional steady-state model of 
the accretion disks produced by WD-NS/BH mergers. Although the 
steady state approximation is clearly invalid just after the merger, 
it represents a reasonable description on timescales ~ r V j sc , during 
which most of the mass accretes at the rate estimated in equation 
(0- We address the stability of our solutions in 94.31 

We adopt a height-integrated model, motivated by previ- 
ous wo rk on RIAFs in th e context of ' advection-dominated' 
models (Narayan & Yil) l994l) with outflows dNaravan & Yilll995l: 
anil 



Blandford & Begelma n 1999). Because RIAFs are geometrically 
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Figure 1. Schematic diagram of accretion and outflows following the 
merger of a WD-NS/BH binary. At the outer edge of the disk (radius 
Rd ~ 10 s 5 - 10 95 cm) the accretion rate M(R d ) is ~ 1(T 4 - 0.1 Af© s~\ de- 
pending on the properties of the merging binary and the viscosity a (eq. (4)). 
The outer composition is He, He-C-O, C-O, or O-Ne, depending on the 
properties of the disrupted WD. As matter accretes to small radii, it ex- 
periences higher temperatures and burns to increasingly heavier elements. 
Outflows are driven from the disk (directly or indirectly) by energy released 
from viscous dissipation and nuclear burning (Fig. 12151 . Elements synthe- 
sized near the midplane are transported to the surface, such that the local 
composition of the wind matches that of the midplane. The outflow com- 
position characterizing each range of radii are shown next to the outgoing 
arrows for the case of a C-O WD (the He WD case is shown in brackets). 
The velocity of the wind is a factor It]^ 2 times the local escape speed, such 
that the inner disk produces the highest velocity outflow. On larger scales 
the ejecta forms a singular homologous outflow with a much lower velocity 
dispersion (see text). At radii r < i?i n ~ 10 6 - 10 7 cm the temperature is suf- 
ficiently high that nuclei are photodissintegrated into 4 He and, ultimately, 
free neutrons and protons. Near the surface of the NS or BH, cooling by 
neutrinos may be important if the accretion rate is sufficiently high. 



thick, averaging over the vertical structure serves only as a crude 
approximation. In standard RIAF models, viscous heating is bal- 
anced by advective c ooling, due either to accretion or outflows (see 
iNaravan et all 1 19981 for a review). As we describe below, a key 
aspect of accretion following WD-NS/BH mergers is the relative 
importance of heating from nuclear reactions compared to viscous 
heating. We term this novel accretion regime a Nuclear Dominated 
Accretion Flow (NuDAF). 

We begin by defining the (positive) steady state accretion rate 



M(r) = -TxrZVr, 



(5) 



where E = 2Hp is the surface density, p is the midplane density, 
H is the vertical scaleheight, and v r < is the radial velocity. We 
allow for the presence of a wind with a mass loss rate M w , which 
we characterize by the local quantity 



p(r) ; 



<91nM w 
dlnr 



(6) 



We do not demand that p(r) be constant with radius, but rather de- 
termine its functional form self-consistently from the solution, as 
described below. 

We assume that the wind exerts no torque on the disk, such 
that it carries away only its own specific angular momentum, 
i.e. j w = ji = r 2 £l, where Q is the angular velocity. This is valid 
provided that the disk is not threaded by a large scale field of suffi- 
cient strength to establish an Alfven radius significantly above the 
disk surface. We furthermore assume that no torque is applied by 
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the inner boundary, such that (at radii much greater than the inner 
disk edge) the radial velocity obeys 

v dQ. 3v 3 (H\ 2 

where we adopt a lShakura & Sunvaevl dl973l) kinematic viscosity 
y = aP/pQ. with a «; 1; H at a /Si < r is the vertical scale 
height; a = (P/p) 1/2 defines the midplane sound speed; and in 
the last expression we assume that Si equals the Keplerian rate 
Si k = (GM/r 3 ) 312 = v k /r to leading order in (H/r) 2 . In adopt- 
ing an anamalous viscosity we have implicitely assumed that an- 
gular momentum by MHD turbulence or gravitational instability 
outwards dominates any inward transport from convection; in 94.41 
we discuss the validity of this point. From equation ((5} we note that 
M = "inVL. 

Mass continuity dM^/dr = dM„/dr implies that 



dlnp dlna 
dlnr dlnr 



P-3, 



while radial momentum conservation 

dv r 



1 dP GM , 
+ - — + — = rQr 



dr p dr 
requires that 

dlnp ^ dlna 
dlnr dlnr 







r 2 (Si 2 - Sll) 



0. 



(8) 



(9) 



(10) 



where we have neglected the first term oc v 2 oc (H/r) 4 . 
The entropy equation reads 



<?iidv 



ds 

-= V ' T Tr 



de dp- 1 
— +P-h— 

or or 



2 Slna 
r \y — 1 dlnr 



dlnp 
dim- 



I dSi \ 2 ^ 9 
[dlnr) ~ 4 



vSi 1 



4 k 



1 + 



V r I \ dlnr dlnr 



of the accreting mass, viz. p 6 [0, 1]. If a fraction of the accreting 
matter does in fact escape, then the upper atmosphere of the disk 
must be preferentially heated. Only in this manner can the wind 
achieve a positive Bernoulli parameter Be w = v 2 /2 > 0, while the 
remaining disk material is accordingly cooled and remains bound 
(Bed < 0), where v w is the asympto tic velocity of the wind. 

Following previous work (e.g. lKohri et al.ll2005l) , we quantify 
the efficiency of wind heating by a parameter tj w = v 2 /2v 2 , which 
equals the ratio of Be w to the local gravitational binding energy. 
This prescription is certainly not unique; for instance, v w could in- 
stead scale with the disk sound speed a. Although for our purposes 
here this distinction is unimportant since and a have similar val- 
ues, the pressure-dependence of wind cooling is different between 
these cases, which has an important effect on the thermal stability 
of the disk ( 943l 

As we describe below, by demanding that the Bernoulli pa- 
rameter in the midplane Bej < and fixing t/ w , this uniquely 
specifies the wind mass outflow rate p(r) (eq. J6)). Although we 
adopt a 'two zone' model (disk + outflow), we do not specify the 
source of wind heating explicitly^ In most of our calculations we 
adopt values for rj„ ~ few, because a terminal speed of the order 
of th e escape speed is a common feature of thermally-driven winds 
(e.g. lLamers & Cassinellil 19991) . 

If the wind is heated, then the disk necessarily cools at the rate 



(11) 



where T is the midplane temperature, y e [4/3,5/3] is the adi- 
abatic index, and .? and e are the specific entropy and internal 
energy, respectively. We have neglected the chemical potential 
term oc yu,-^ because all species (including electrons) are non- 
degenerate at the radii where nuclear burning commences. Equa- 
tion dl It may be interpreted as the balance between "cooling" 
from advection through the disk g a[lv and others sources of heat- 
ing/cooling q = q visc + q nuc + q v , including those resulting from 
nuclear reactions q nuc (p, T) (see eq. 1141 below), winds q w (see be- 
low), and viscous dissipation 



,(12) 



where we have used equation 010) in the last expression and neglect 
terms oc d(H/r)/dr. We again note that radiative cooling may be 
neglected because the timescale for photon diffusion from the disk 
midplane is much longer than the viscous time time (eq. (3))- In 
94.41 we describe a scenario in which this might be violated, but 
this alternative picture of the disk evolution is radically different 
from the model developed in this section. 

An important characteristic of RIAFs without outflows (p = 
0) is that the Bernoulli parameter 

1,1,, P GM 

Be d = -v 2 + -r 2 Q. 2 + e+ (13) 

2 2 p r 

is generically positive dNaravan & Yilll994l) . This implies that, in 
principle, material in the midplane has sufficient energy to es- 
cape to infinity. This fact has been used pre viously to argue that 
RIAFs are susceptible to pow erful outflows dNaravan & Yilfl995l ; 
iBlandford & Begelman|[l999l) that carry away a substantial portion 



(yl/2 - Be d ) dM A 



lixLr 



31nr 



--p(n„ -Be>n 2 , 



(14) 



where Be d = Be,j/v? is the 'normalized' disk Bernoulli function. 
By combining equations d 1 It— dT4t we find that 



(!)' 



3y — 1 <91na 1 <91np 



y - 1 dlnr 2 dlnr 



2 ^/niic 

"31*2; 



•Be d )-- (15) 



Xaif 



(16) 



Finally, nuclear reactions change the composition of the ac- 
creting material according to 

dX A 
V ' — 

where Xa is the mass fraction of isotope with mass number A, and 
X A \ P is the nuclear reaction rate. Equation 1 1161 shows that material 
burns for approximately the local accretion timescale ~ r/v r at any 
radius. We assume that burning occurs at constant pressure because 
in general the burning timescale is long compared to the dynamical 
timescale over which vertical pressure balance is established. 

We neglect the effects of convective/turbulent mixing in the 
disk, which may act to smooth radial abundance gradients; mix- 
ing could, in principle, be modeled by including an additional term 
K v m [ x V 2 X A to the right hand side of equation d 1 6b . where is 
the diffusion coefficient. Neglecting mixing is justified as a first ap- 
proximation because compositional changes typically occur over 
a radial distance > H. Studies of the diffusion of passive contami- 
nants in numerical simulations of the MRI furthermore demonstrate 
that < v (e.g. lCarballido et"al . 2005), such that it appears un- 
likely that burned material will diffuse upstream. 

The nuclear reaction rates X A (p, T) and nuclear heating rate 
<7nuc(P> T) in equations 1 11 It and j Kit are calculated using a 19 iso- 
tope reaction network that includes ff-capture, heavy-nuclei, and 

' One possible source of coronal heating is the dissipation of MHD waves. 
Waves may be excited by turbulence due to the MRI or vertical convection, 
the latter of which is likely due to the strong temperature dependence of the 
nuclear heating rate. 

2 See http : //www . cococubed . com/code_pages/burn_helium . shtml 
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(a,p)(p,a) reactions dTimmesll999h . The temperature T(r) and adi- 
abatic index y(r) are calculated using a standard equation of state, 
which includes ideal gas, radiation, and degeneracy pressure. 

We calculate solutions by integrating equations l[8}, dl5t , 
from the outer boundary r = R4 inwards, in order to obtain a(r), 
p(r) and X A (r). The wind mass loss rate p(r) (eq. (6j) is determined 
by demanding that the Bernoulli function is regulated by outflows 
to a fixed value Be d = -0.2 - (marginally bound disk). Boundary 
conditions are determined by fixing the outer accretion rate M d (^ d ) 
(eq. I?)) and by assuming that the inflowing material reaches a self- 
similar evolution with a oc Vk oc r~ 1/2 , such that ^|« d = -1/20 We 
terminate our solutions at the radius R iu < 10 7 cm at which heavy 
nuclei are photo-disintegrated into free nuclei. Interior to this point, 
neutrinos from the reactions e~ + p — > n + v e and e + + n — > p + v ? 
become an important sourc e of cooling (e.g. iNaravan et al.ll200ll ; 
IChen & Belobor odovll2007h . Depending on M and a, the disk may 
become radiatively efficient and geometrically thin, such that strong 
outflows may be suppressed close to the surface of the central 
NS/BH. Figure [T] shows a schematic diagram of the structure of 
the disk and outflows. 

The composition at the outer boundary r = R^ is that of the 
tidally disrupted WD. For C-O WD s we adopt an initi al compo- 
sition X l2 = 0.5 and X t6 = 0.5 dSalaris etal .111 9971) . We also 
pe rform calculations for O-Ne WDs (X l6 = 0.7; X 20 = 0.3; 
e.g. iGutierrez et al.ll2005l) hybri d He-C-0 WDs (X 4 = 0.2; X u = 
0A-.X lf , = 0.4; Irian etalJboOCl) , and pure He WDs (Xt = 1), al- 
though the latter may occur less frequently in Nature due to the 
limited parameter-space available for unstable mass transfer (Sj2j. 
Although we focus on WD-NS binaries, we also consider 3M Q BH 
companions. Table Q] enumerates the models calculated in the next 
section. 



C-0 WD/NS, M(R d ) = 3x1CT 3 M Q s , tj w =2, a=0.1 




(a) Radial Profile of Disk Properties 




(b) Heating/Cooling Terms in the Entropy Equation 



4 RESULTS 

Figures l2T5l and TablesQ~]and[2]summarize the results of our calcu- 
lations. In this section we describe the disk structure Q4.lt . outflow 
composition (j j4,2> , and thermal stability of our solutions ( 94.3t . In 
94.41 we discuss the possible effects of convection on our results 
and in 94.51 we present an analytic estimate of the importance of 
nuclear burning. 



4.1 Disk Structure 

Figure [2] shows our baseline model for accretion following the 
merger of a 0.6M o C-O WD with a 1.4M S NS. In this exam- 
ple we adopt characteristic values for the initial accretion rate 
M(Ri) = 3 x 1O~ 3 M0 s -1 , viscosity a = 0.1, disk Bernoulli pa- 
rameter Be d = -0.1 (eq. 11131 ), and wind strength r] w = 2. The top 
panel shows the radial profiles of quantities characterizing the disk 
thermodynamics. Note the following: (1) the disk is geometrically 
thick, with aspect ratio H/r w 0.4 at all radii, as expected for RI- 
AFs; (2) gas pressure P ga8 exceeds radiation pressure P m a, although 
they are comparable at small radii < 10 7 cm; (3) M decreases by a 
factor -10 between R a as 10' cm and the inner radius R m x 10 7 
cm, such that > 80% of the accreting mass escapes in outflows. This 



Regardless of the precise outer boundary condition adopted, the self- 
similar solution rapidly obtains jNaravan & Yi 1994). Self-similarity is 
maintained until nuclear burning becomes important, which introduces an 
additional scale into the problem. 
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(c) Composition Profile 

Figure 2. Accretion following the merger of a 0.6M o WD with a 1.4M Q 
NS. We adopt fiducial values T] w = 2 for the wind efficiency, Be d = -0. 1 
for the Bernoulli function of the disk (eq. 1131 ). and a = 0. 1 for the disk 
viscosity. The latter translates into an accretion rate M(/? d ) =a 3 X 10~ 3 M o 
s~' at the outer edge of the disk = 10 9 cm (eq. (4))- Top: Disk tem- 
perature T (Blue), surface density E (Red), aspect ratio H/r (Turquoise), 
ratio of radiation to gas pressure P ra d/Pgas (Green), negative of the electron 
chemical potential yi e divided by kT (Pink), accretion rate M normalized to 
M(Rd) (Brown), and wind mass loss index p (Orange; eq. (6)), as a func- 
tion of radius. Middle: Sources of heating and cooling (eq. 1111 ). in ratio to 
the viscous heating rate q vlsc : nuclear reactions ij nuc (Solid Blue); advective 
cooling within the disk q.^ = Tv r (ds/dr) (Solid Dotted); and wind cooling 
q w (Dashed Green; eq. 1141 ). Bottom: Mass fraction Xa of various isotopes 
in the disk midplane. 
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He WD/NS, M(R d )=1(T 3 M s , tj w =2, a=0.1 



He-C-0 WD/NS, M(R d )=1CT 3 M s~ 1 , t? w =2, a=0.1 
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(b) Heating/Cooling Terms in the Entropy Equation 
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(c) Composition Profile 

Figure 3. Similar to Figure [5] but calculated for a pure He WD. Values of 
r] w = 2, Be d = —0.1, and a = 0.1 are the same as the C-0 case, but the 
latter corresponds to a lower accretion rate M(Ri) 10~ 3 M Q s - ' due to the 
larger outer radius i?d = 2 X 10 9 cm. 



behavior is reflected in the wind mass loss index p (eq. f6]), which 
fluctuates with radius, depending on the local heating/cooling from 
nuclear reactions. 

The middle panel in Figure [2] shows sources of heating and 
cooling that contribute to the entropy equation (eq. Ill 11 ). in ratio to 
the viscous heating rate q vlsc (eq. 1121 ). At large radii (r > 5 x 10 7 




(a) Radial Profile of Disk Properties 




(b) Heating/Cooling Terms in the Entropy Equation 




(c) Composition Profile 
Figure 4. Same as Figure[5] but calculated for a 'hybrid' He-C-O WD. 



cm), both viscous dissipation and nuclear burning contribute to the 
heating, while cooling results from advection to smaller radii </ adv 
and mass loss due to winds q w (eq. 11141 ). The nuclear heating rate 
achieves local maxima at radii where the temperature reaches the 
threshold necessary to burn the next element (e.g. carbon burns at 
r at 3 x 10 8 cm and T ~ 1.5 x 10 9 K). At these locations the energy 
generated by nuclear reactions is significant compared to that re- 
leased by viscous dissipation; this general behavior has important 
implications for the thermal stability of the disk ( 94.3t . At smaller 
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O-Ne WD/BH, M(R d ) = 3x1CT 2 M s~\ 7^=2, a=0.1 
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(a) Radial Profile of Disk Properties 
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(b) Heating/Cooling Terms in the Entropy Equation 
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(c) Composition Profile 

Figure 5. Similar as Figures I2l4l but calculated for the merger of a O-Ne 
WD with a 3M Q BH. We adopt values of 7] w = 2, Be' d = -0.2, and a = 0.1, 
and M(R A ) = 3 X 10" 3 M Q s _I , R d = 5 X 10 s cm. 



radii (r < 5 x 10 7 cm), by contrast, q nuc < instead contributes 
to cooling the disk. This occurs once temperatures are sufficiently 
high (T > 6 x 10 9 K) to photo-disintegrate nuclei endothermically . 

The bottom panel of Figure [2] shows the mass fractions of in- 
dividual elements X A as a function of radius. The composition at 
the outer edge of the disk is that of the initial WD: half 12 C and half 
16 0. Moving to smaller radii, carbon first bums to 20 Ne and 4 He, the 



latter of which rapidly capture onto 16 0, eventually forming 24 Mg. 
At smaller radii and higher temperatures, 16 burns to 28 Si. Mov- 
ing yet further in, photo-dissociation releases a-particles, which 
immediately capture onto heavier nuclei ('Si burning'), producing 
32 S, 36 Ar, 40 Ca, 56 Ni, and 54 Fe. Finally, at the smallest radii, 54 Fe 
photo-disintegrates into 4 He, which itself then disintegrates into 
free protons and neutrons. 

In Figure [3] we show similar results to Figure [2] but for the 
case of a 0.3M o pure He WD. We adopt the same values for t] w = 2, 
Be d = -0.1, and a = 0.1, but the latter corresponds to a lower ac- 
cretion rate M(R d ) ^ 10 _3 M o s _1 than in the C-0 WD case, due 
to the larger outer radius /Jj = 2 x 10 9 cm. The qualitative disk 
structure is similar to the C-0 case, but a few differences should 
be noted. First, mass loss is considerably greater from the outer 
portions of the disk (p > 1). This is due to the substantial energy 
released by He burning and its low threshold temperature; burn- 
ing thus occurs at large radii where the gravitational potential is 
shallow. Mass loss greatly reduces the midplane density, which in- 
creases the importance of radiation pressure. Because the tempera- 
ture is lower in radiation-dominated than in gas pressure-dominated 
disks, burning heavier elements is delayed until smaller radii than 
in the C-0 case. Also note that the production of 16 0, 20 Ne and 
24 Mg are largely bypassed (e.g. in comparison to the C-0 case) 
because additional a captures are rapid following the onset of the 
[rate-limiting] triple-cf reaction. 

Figure [4] shows results for a 0.4M o 'hybrid' He-C-O WD, 
again for an accretion rate M(Ra) = 10~ 3 M o s -1 . Although the inner 
structure of the disk resembles the C-0 case (Fig. [2}, an important 
difference is again the onset of He burning at relatively large radii 
(low temperature), this time via the reaction 4 He+ 16 0— > 20 Ne+y. 
As we discuss in 94.31 the sensitive temperature dependence of this 
reaction may render such disks thermally unstable. 

Figure [5] shows results for a 1.2M Q O-Ne WD accreting onto 
a 3M Q BH. The compact outer radius of the disk in this case R^ =; 
5 x 10 s cm results in a high accretion rate M{R A ) - 3 x 10~ 2 M o s _1 
for a = 0.1. Without He or C present, the first element to bum in 
this case is 20 Ne at T « 2 x 10 9 K, but otherwise the composition 
profile is similar to the C-0 case (Fig. [2}- 

In addition to the calculations presented above, we have ex- 
plored the sensitivity of our results to variations in the binary pa- 
rameters, WD composition, viscosity a, and wind mass loss pa- 
rameter rj w (see Table [2]l . In general, we find that increasing t] w at 
fixed Be d decreases the amount of mass loss because wind cooling 
is more effective. By contrast, changing a (along with a compen- 
sating change in M(R&)\ eq. |4)) has a smaller effect. Our results are 
also robust to realistic variations in the composition of the disrupted 
WD. 

Notably absent from the discussion thus far is the merger of a 
massive C-0 or O-Ne WD with a NS companion. The large mass 
ratio q ~ 1 in this case results in a compact disk (R^ ~ 10 8 cm) with 
a very high initial accretion rate M(Ra) ~ 0.1 M G s _I for a = 0.1. 
The initial temperature of the disk is sufficiently high that nuclear 
burning begins already at r ~ Ra during the circularization pro- 
cess itself. Since virialization has not yet occurred, the WD mate- 
rial forming the disk may still be degenerate, potentially leading to 
an explosive situation. Furthermore, as we discuss in 94.31 accre- 
tion in these systems may be thermally unstable, even under non- 
degenerate conditions. Clearly, many of the simplifications adopted 
in our model break down in the case of massive WD-NS mergers. 
Obtaining even a qualitative understanding of such events may re- 
quire multi-dimensional simulations that include the effects of nu- 
clear burning, even during the merger process itself. 
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4.2 Outflow Properties 

A schematic diagram of the outflows from WD-NS/BH merger 
disks is shown in Figure [T] The total mass loss rate from the 
disk is approximately given by M w =; M{R A ) - M(R m ), where 
R m ~ 10 6 5 - 10 7 5 cm is the radius interior to which the disk is 
primarily free nucleons. We adopt this inner boundary because if 
neutrino cooling is important at radii r < R m , then additional mass 
loss will be suppressed (although in j )5.2| we discuss the possibility 
that the ejecta does contain a small fraction of free nucleons). If 
most of the WD accretes at the characteristic rate M(Ra) (eq. l4l). 
then the total ejecta mass is M e j ss / c jM WD , where / c j = J is 
the fraction of the accreted mass lost in winds. Depending on t/ w , 
we find that / w ~ 0.5 - 0.99[0.90 - 0.99] for C-0[He] WDs (Table 
0, corresponding to M cj ~ 0.3 - 1.3[0.3]M o . 

The radial structure of the disk is imprinted on the composi- 
tion of the outflows (Fig.[TJ. Table [2]provides the mass fractions in 
each element of the wind ejecta, which are calculated as 



-f 

M w Jr 



Rd dhl 
X-A—z—dr 
or 



(17) 



Here we have assumed that the local surface composition of the 
disk is similar to that in the midplane due e.g. to turbulent mix- 
ing. To compare with our C-O disk models, we also list the 
abundances predicted by the well-studied 'W7' model for SN la 
from the delayed-de tonation of a Chandrasekhar-mass C-O WD 
dNomotoetaflll984h . To compare with our He disk models, we 
list the range in ejecta abundances predicted from models for un- 
st able He shell bur ning on the surface of a C-O WD, as calculated 
bv lShen et aljfeoid) for different values of the mass of the WD and 
He layer. 

In our fiducial C/O model (Fig.O the majority of the ejecta 
is unburnt 12 C (26%) and 16 (39%), with the remainder in 28 Si 
(10%), 24 Mg (7%), 4 He (6%), 20 Ne (5%), 54 Fe (3%), 32 S(2%), 56 Ni 
(1.1%), 36 Ar (0.4%) and 40 Ca (0.3%). Compared to the W7 SN la 
model, the fractional abundances of unburnt C/O, 20 Ne, and 24 Mg 
are significantly higher. Despite these differences, the abundances 
of the elements Si, S, and Ca, generally responsible for producing 
the most prominent spectral features in normal la, are similar to the 
W7 model within a factor < 2-3. Importantly, the yield of radioac- 
tive 56 Ni is substantially lower than in normal la SN ejecta; this has 
important implications for the luminosities of optical transients as- 
sociated with WD-NS/BH mergers ($52\. 

In our fiducial He model (Fig. [3} the vast majority of the 
ejecta is unburnt 4 He (63%), with the remainder in 12 C (15%), 
28 Si (10%), 32 S(8%), 36 Ar (2%), 54 Fe (0.8%), 24 Mg (0.6%), 40 Ca 
(0.3%), 20 Ne (0.2%), 16 (0.2%), and 56 Ni (0.1%). Although the 
unburnt 4 He fraction is similar to the He shell burning models of 
Shen et al. (2010), the predicted abundances of the radioactive iso- 
topes 44 Ti, 48 Cr, and 56 Ni are all significantly lower, while those 
of intermediate mass elements are significantly higher. In addition 
to the above elements, disk outflows may in all cases also contain 
a tiny fraction of 'H from the very inner disk, the implications of 
which are briefly discussed in 35.21 

Winds from smaller radii in the disk achieve higher velocities 
V w oc vt k r~ </2 than those from further out. Nevertheless, most 
of the ejecta will ultimately reside in a single 'shell' of material 
with a much lower velocity dispersion. This is because, initially 
after the disruption (on timescales t <K ? v j sc ; eq. (3)), the disk is 
concentrated at large radii = R A . Slow outflows from this early stage 
will thus have time to encase the system before outflows have even 
begun in earnest from the inner disk on a timescale t > f v j sc . Fast 



ejecta from the inner disk will thus collide with the slower material 
ejected prior. Because this interaction conserves energy (radiative 
losses are negligible at early times), the ejecta will achieve a mean 
velocity v c j, which may be estimated by averaging over the total 
kinetic energy of winds from each disk annulus: 



1-2 1 f Rj <3M W 2 , 



(18) 



Table [2] shows that for C-O WD/NS mergers v w is typically ~ 
1 - 3 x 10 4 km s -1 , while for He[0-Ne] mergers v w is somewhat 
lower[higher]. 

The reason that v w can exceed the typical velocity of SN 
la ejecta is the deeper gravitational potential well in the case of 
NS/BH accretion. As a toy example, if mass loss were to occur 
with a constant power-law index p = const < 1 (eq. f6]), then 
integrating equation (18t gives 



[ l-p R in \ Si I J ' 

h"f'*(e 1/2 > 



p<\ 
p=l 



(19) 



where in the first expression we have assumed that R A » R^. A 
larger wind efficiency factor t/ w thus increases v c j, both because 
v w oc rjlJ 2 and because larger r/ w decreases the mass loss rate p(r) 
required to cool the disk. Equation d 1 9 1 shows that (for t]„ = 1) the 
minimum value of v c j (p ~ 1) is of the order of the orbital veloc- 
ity at r ~ R d . This varies from « 3 x 10 3 km S" 1 for He WD-NS 
mergers to « 10 4 km s _I for O-Ne WD-NS/BH mergers. 



4.3 Thermal Stability 

The calculations presented in the previous sections assumed steady- 
state accretion. This is not valid, however, if the disk is thermally 
unstable, which might be expected due to the sensitive temperature 
dependencies of nu clear reactio ns. The criterion for thermal stabil- 
ity may be written dPirarj[l978l) 



dlnq + 



dlnH 



dlnq 



dlnH 



d\nq + 



dlnP 



dlnq 



dlnP 



(20) 



where q + = q v i x + q„u C and if = q^ + q w are the total heating 
and cooling rates, respectively (see eq. Ill 11 ), and the second con- 
dition follows because P = pa 2 oc Q^LH. Viscous heating obeys 
<?visc K R 2 , while advective cooling obeys g a( j v oc q vlsc (Hlr) 2 oc P 4 . 
In our model, we have assumed that the terminal speed of disk 
winds v w is proportional to the local escape speed oc y k (eq. 1141 ). 
such that q w oc g visc oc P 2 . However, the precise mechanism respon- 
sible for driving the wind is uncertain; we could equally well have 
assumed that v w is proportional to the sound speed a, in which case 
q w oc (///r) 2 V| c oc P 4 . Considering both cases, the stability criterion 
in equation j20\ can equivalently be written as an upper limit on 
the fraction / nuc = q nuc /q + of the disk heating supplied by nuclear 
reactions, viz. 



,/nuc < - 



2»-/») _ f 
„_2 — Jniic.thl » 

/Iiiic,th2, 



2 
n-2 



Vw K V esc 
V w oc a, 



(21) 

where n = f a ^^ c J | E and f w = q„/q~~ is the fraction of the cooling 
resulting from winds. Equation (121) shows that accretion is stable 
at a given radius if either (1) nuclear burning is absent entirely; or 
(2) the burning rate has a weak pressure dependence n < 4 - 2/ w or 
n < 4, in the cases v w oc y k and v w cc a, respectively. 
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Table 1. Properties of Disk Models 



model 


(Mq) 


(M Q ) 


(cm) 


M(R A ) 
(M Q s" 1 ) 


a 


Be/ 5 

d 






M ■ , W) 
M(«d) 


(10 4 kms _1 ) 


unstable?"' 1 


NS_C-0_1 


1.4 


0.6 


10'' 


3 X 10 3 


0. 1 


-0. 1 


2 


X\2 = 0.5, X\ 6 = 0.5 


0.84 


2.8 


no 


NS_C-0_2 


1.4 


0.6 


10 l) 


3 x 10~ 4 


0.01 


-0.1 


2 


- 


0.86 


2.8 


no 


NS_C-0_3 


1.4 


0.6 


10' 


3 x 10~ 3 


0. 1 


-0.2 


1 




0.99 


1 .6 


no 


NS_C-0_4 


1.4 


0.6 


10 


3 x 10 -3 


0. 1 


0.0 


3 




0.57 


3.2 


no 


NS_C-0-He_l 


1.4 


0.4 


2 x 10 J 


10 - 


0. 1 


-0.1 


2 


X if> = 0.4, Xn = 0.4, Xi = 0.2 


0.95 


2.5 


yes 


NS_C-0-He_2 


1.4 


0.4 


2 x 10 9 


10" 4 


0.01 


-0.1 


2 




0.95 


2.4 


yes 


NS.O-Ne (/> 


1.4 


1.2 


2x 10 8 


io-' 


0.1 






X 16 = 0.6, X 20 = 0.4 






yes 


BH.C-OJ 


3 


1.0 


5 x 10 s 


3 x IO" 2 


0.1 


-0.2 


2 


X n = 0.5, X m = 0.5 


0.72 


3.6 


yes 


BH.C-OJ 


3 


1.0 


5 x 10 8 


3 x 10~ 3 


0.01 


-0.2 


2 


check 


0.72 


3.6 


yes 


BH_C-0_3 


3 


1.0 


5 x 10 s 


3 x IO" 2 


0.1 


-0.25 


1 




0.96 


2.3 


yes 


BH.O-Ne.l 


3 


1.2 


5 x 10 8 


3 x 10~ 2 


0.1 


-0.2 


2 


X u = 0.6, X 10 = 0.4 


0.83 


4.6 


yes 


BH.O-Ne.2 


3 


1.2 


5 x 10 8 


3 x 10~ 3 


0.01 


-0.2 


2 




0.72 


3.5 


yes 


NS-He.l 


1.2 


0.3 


2x 10 9 


io- 3 


0.1 


-0.1 


2 


X 4 = 1 


0.95 


2.2 


no 


NS-He_2 


1.2 


0.3 


2x 10 9 


io- 4 


0.01 


-0.1 


2 




0.96 


2.0 


no 


NSJie.3 


1.2 


0.3 


2x 10 9 


io- 3 


0.1 


-0.2 


1 




0.99 


1.0 


no 



{a) Mass of the central NS or BH. ''''Outer edge of the disk R d = SrlofO + <?)~ 3 , where q = My/o/M and Rrlof is the binary separation at RLOF (eq. 0). 
^'Normalized Bernoulli parameter in the disk (eq . 11131 ). Ratio of the specific kinetic energy of the wind from any radius to the local gravitational binding 
energy. ' e 'lnitial composition of the disk from the disrupted WD. '^'Ratio of the total mass loss rate in the wind M w = M(R L \) - M(R m ) to the initial accretion 
rate M(Ri) at the outer edge of the disk. ' s 'Mean velocity of the ejecta from disk winds (eq. 1181 ). ''''Whether the disk is thermally unstable due to nuclear 
burning at any radius, according to the most stringent criterion / nuc < / nu c,th2 given in equation j20\ . '''Characteristic values shown for illustration only. The 
disk model presented in this paper cannot be applied to the merger of massive WD-NS binaries because nuclear burning begins already at r » Rj (see j|4.2E . 



Table 2. Elemental Mass Fractions in the Wind Ejecta X w a (eq. 1171 ) 



mode] 


4 Hc 


12 C 


L6 Q 


20 Ne 


24 Mg 


28 Si 


32 S 


36 Ar 


4l) Ca 


44.fi 


48 Cr 


52 Fc 


54 Fe 


56 Nj 


W7 Ia<°> 


0.0 


0.036 


041) 


0.0014 


0.0094 


041 


0.060 


0.1 1 


0.0086 


5.7 X Ur'< 


3.5 X 10"' 




0.054 


0.42 


NS-C-O-1 


0.063 


0.26 


0.39 


0.049 


0.065 


040 


0.023 


0.0043 


0.0031 


2.1 x ur 5 


1.0 X IO -4 


9.7 X 10" 4 


0.032 


0.011 


NS_C-0_2 


0.052 


0.26 


0.38 


0.058 


0.067 


0.096 


0.021 


0.0036 


0.0025 


1.3 X 10 -5 


7.0 X I0 -5 


8.0 X 10" 4 


0.042 


0.011 


NS_C-0_3 


0.0091 


041 


0.47 


0.029 


0.034 


0.034 


0.0068 


0.0012 


8.0 x ur* 


3.7 x ur 6 


1.4 X 10" 5 


1.0 X IO -4 


0.0061 


1.0 X 10-' 


NS.C-O-4 


0.070 


0.21 


0.36 


0.057 


0.076 


0.13 


0.029 


0.0055 


0.0040 


2.6 X 10~ 5 


1.3 x Ur 4 


0.0012 


0.035 


0.014 


NS_C-0-He_l 


0.12 


0.33 


0.18 


0.093 


0.16 


0.087 


0.0064 


0.0010 


6.5 x ur 4 


4.2 X 10-" 


1.6 x Ur 5 


1.2 x IO" 4 


0.010 


0.0010 


NS_C-0-HeJ 


0.1 I 


0.33 


0.17 


0.095 


0.19 


0.071 


0.0054 


9.5X10- 4 


7.1 x ur 4 


4.5 X 10-* 


2.3 X 10-' 


2.3 X IO -4 


0.011 


0.0027 


BH.C-O.1 


0.15 


6.7 X ur 4 


0.20 


0.043 


0.11 


0.24 


0.047 


0.0081 


0.0063 


4.8 x ur 5 


2.9xH)- 4 


0.0036 


0.12 


0.053 


BH.C-O.2 


(1.14 


5.7 X IO -5 


04 9 


0.021 


0.098 


0.25 


0.047 


0.0062 


0.0039 


1.2 x ur 5 


6.3 X I0 -5 


7.4 x 10" 4 


0.22 


0.011 


BH_C-0_3 


0.062 


0.12 


0.30 


040 


0.098 


047 


0.032 


0.0053 


0.0040 


3.0 X 10 -5 


1.8 X IO -4 


0.0021 


0.065 


0.031 


BH_ONc_l 


0.15 


4.7 x 10-" 


0.30 


0.045 


0.035 


047 


0.076 


0.016 


0.013 


5.5 x ur 5 


3.0 X IO -4 


0.0036 


0.12 


0.054 


BH_0-Ne_2 


0.13 


2.8 X IO -6 


0.27 


0.024 


0.037 


04 8 


0.080 


0.015 


0.011 


1.8 X 10~ 5 


6.6 x Ur 5 


7.1 x ur 4 


0.22 


0.010 


ShenlO .la" 1 


0.30 - 0.45 














6 - 700 X 10" 5 


0.001 -0.04 


0.007 - 0.1 


0.01 -0.1 


0.03-0.1 




0.1-0.7 


NS_He_l 


0.63 


04 5 


0.0022 


0.0017 


0.0055 


0.095 


0.082 


0.018 


0.0027 


7.1 X 10-* 


1.8 X 10" 5 


l.i x ur 4 


0.0082 


9.7 x ur 4 


NS_He_2 


0.60 


0.23 


0.0031 
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Figure [6]shows / nllc (solid line), / nuc>1 hi (solid lines), and/ nuc th2 
as a function of radius, calculated for our solution from Sj4] corre- 
sponding to accretion following a 0.6M o C-0 WD (Fig. [2j- Fig- 
ure [6] shows that in the case v w ec a, the disk is stable at all radii, 
i-e. /nuc < /nuc,th2- If, however, we instead assume that v w cc y k , then 
we find that the disk is unstable (f mc > f ttll ) at the radii ~ few x 10 8 
cm where carbon burning peaks. Although we do not show this case 
explicitly, we find a simlilar result for He WDs: the triple-a reac- 
tion is not sufficiently pressure-sensitive for instability if v w cc a 
(i.e. n < 4), but for v w cc v k the disk is unstable at radii ~ 3 - 4 x 10 s 
cm where He burning peaks. We conclude that whether C-O or He 
WD disks are thermally stable depends on the pressure dependence 
of wind cooling, which in turn depends on the [uncertain] mecha- 
nism responsible for driving the wind. 

In the case of O-Ne and hybrid WDs, thermal instability seems 
more assured. In the bottom panel of Figure [6] we show / luc and 



/nuc,th2 ( e q- 1211 ) for cases corresponding to the merger of a 0.4M o 
He-C-0 WD with a NS (Fig.©; a 1.2M G O-Ne WD with a 3M Q 
BH (Fig. [5)1; and a 1.0M o C-O WD with a 3M G BH. In all three 
cases we find that over some range of radii, / nuc exceeds even 
the more conservative threshold for instability / nuc _ t h2- This implies 
that our calculated solutions are thermally unstable, indpendent of 
whether v w cc v k or v w cc a. Although addressing the implications 
of unstable burning is beyond the scope of this paper, we speculate 
that the end result may be a complicated, time-dependent evolu- 
tion. One possibility is 'limit-cycle' behavior, in which the initial 
inflow is halted by runaway burning, before resuming again later 
with fresh fuel (a possi ble analog are d warf nova outbursts of cat- 
aclysmic variables; e.g. lCannizz 5 fl993l Though it is not obvious, 
the steady-state solutions constructed here may remain a reasonable 
description of the average flow over many cycles. 
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Figure 6. Analysis of the thermal stability of the accretion solutions pre- 
sented in j|4. II Quantities shown are the fraction of the total heating due to 
nuclear burning f lmc = q nuc /q + (dotted line) and the threshold fractions for 
thermal instability / nU c,thl (dashed line) and / nuc th2 (solid line) under the 
assumption that v w oc v* and v w oc a, respectively (eq. 1211 ). Top : Solu- 
tions corresponding to accretion following the merger of a 0.6M o 1.4M e 
WD with a NS at the rate M(R d ) = 3 X 1(T 3 M G s" 1 (Fig.0. Note that al- 
though f nuc < f aucAh2 at all radii, / nuc exceeds / nuc , thI across some range in 
radius. This shows that whether the disk is in fact stable depends on the de- 
tailed pressure-dependence of wind cooling. Bottom: Accretion following 
the merger of a 1M G C-O WD with a 3M Q BH (Black); 1.2M Q O-Ne WD 
with a 3 M© BH (Blue); 'hybrid' He-C-O WD with a 1 AM e NS (Red). In all 
three cases f mc > / nuc t h2 over some range in radii. The solutions are thus 
unstable to runaway heating from nuclear burning, regardless of whether 
v w oc vt or v w oc a. 



4.4 Effects of Convection 

It is well known that ADAFs a re unstable to rad ial con- 
vection jNaravan& Yil 1 19941 : llgumenshchev et al.l l2000l : 
lOuataert & Gruzinovl l200fj|) . In the case of NuDAFs, the disk 
may also be unstable to vertical convection because the nuclear 
heating peaks sharply near the midplane due to its sensitive 
temperature dependence. Although our calculations have thus 
far neglected its effects on the transport of angular momentum, 
energy, and composition, in this section we describe how strong 
convection could alter our conclusions. 

If convection transports angular momentum outwards -similar 



to MHD turbulence produced by the MRI-then it would act simply 
to enhance the effective value of "or", which we have shown does 
not qualitatively affect our results. If, on t he other hand, convectio n 
transports angular momentum inwards dRvu& Goodman! Il992l) . 
then an qualitatively different type of solutio n could obtain: a "con- 
vection dominated accretio n flow" (CDAF; lOuataert & Gruzinovl 
200G; [Naravan et all 2000). In CDAFs, the angular momentum 
transported outwards by viscosity is balanced by the inward trans- 
port by convection, such that the net mass accretion rate is very 
small. Because in steady state the convective energy flux F c oc 
pc^r 2 is constant with radius, CDAFs are characterized by a more 
shallow density profile p oc r~ I/2 than for normal ADAF solutions 
(p oc r~ 3/2 ). If present, the CDAF would extend from the inner edge 
of the RIAF at /{;„ ~ 10 6 - 10 7 cm to much larger radii, where the 
outflowing energy is released through radiation or outflows. 

Coincidentally, the density profile of a CDAF scales the same 
way with radius as in ADIOS (wind) models with "maximal" mass 
loss (p = 1; M oc r). Our steady-state calculations with high mass 
loss rates (low rj w ) thus may also be applied to the CDAF case. 
Note, however, that the physical inteipretation is completely dif- 
ferent: in wind models, the outflowing mass escapes the system 
entirely, while in CDAFs it remains bound and is simply 'recycled' 
by mixing to large radii, from which it accretes again later. 

If the outer edge of the disk remains near the initial radius 
of the disrupted WD Rj, then convection simply increases the 
timescale for matter to accrete from that given in equation (0 by 
a factor ~ ~ ^ ~ 10-1 3 , or from minutes to hours or 

days. The radial structure of the disk is similar to the maximal wind 
case, but the observational signature of the event as described in §5\ 
will be altered. With no outflows present, the central NS or BH ac- 
cretes significantly more mass, potentially resulting in more ener- 
getic high energy emission ( 35.11 ). On the other hand, with little or 
no radioactive material ejected, the SN-like optical counterpart (if 
any) would be much dimmer (j j5. 2b . Outflows may still occur from 
the outer radius of the CDAF in this scenario, but the mean velocity 
is likely be much lower and the 56 Ni fraction would be smaller. 

A more radical possibility is that the outward convective en- 
ergy flux will 'feed back' on the outer disk, alterating its dynam- 
ics entirely. If pressure forces cause the outer reservoir of mass to 
expand significantly from its initial radius ~ R A , then rotational 
support will no longer be important. As a much larger, quasi- 
spherical accreting envelope, the structu re may thus come to re- 
semble a similar 'T horne-Zytkow object' dThorne &Zvtkowll9'75l) 
or 'quasi-star' (e.g. Be gelman et al. 2008 1. In this case, the accre- 
tion timescale will be even long than estimated above, because it 
instead depends on the rate that the convective energy escapes the 
outer boundary in radiation or winds. A quantitative exploration of 
such a model is an important and interesting excercise for future 
work, but is beyond the scope of this paper. 

Despite the important possible effects of convection discussed 
above, it is not clear that CDAFs are in fact relevant to WD-NS/BH 
mergers. As already noted, whether the flow is best described as a 
CDAF or ADAF (with winds) depends on the strength (and direc- 
tion) of convective angular momentum transport relative to other 
torques. Although radial convection in hydrodynamic simulations 
of RIAFs indeed appears to balance outwards viscous transport for 
a < 0.05 (Naryan et al. 2000; Narayan et al. 2004), higher val- 
ues of a ~ 0.1 are generally inferred from o bservations of ion ized 
disks in a variety of astrophysical contexts dKing et al.ll2007h . In 
WD-NS/BH mergers, angular momentum transport due to gravita- 
tional instabilities may also compete with convection. Finally, it has 
been argued that global MHD simulations of RIAFs more closely 
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resemble the ADIOS picture of outflows em ployed in this paper 
dHawlev et alll200ll ; lHawlev & Balbusll2002h . although this con- 
clusion may depend on the [uncertain] saturation strength of the 
MRI and the ultimate fate of the mass leaving the grid. 



4.5 Analytic Estimate of the Importance of Nuclear Burning 

We conclude this section with an analytic estimate of the impor- 
tance of nuclear burning in hyper-accreting disks, i.e. we address 
when a disk is in the NuDAF regime. Focusing on the burning 
of one isotope A, the nuclear heating rate may be written q nuc = 
Aenuc/tfJnfbu,,,, where Ae nuc = XaQ/A is the nuclear energy released 
per nucleon, Q is the Q-value of the reaction, m„ is the nucleon 
mass, f bum = Ar bum /v r (/burn) is the time spent burning at radius 
t = rbum, v r = 3v/2r is the accretion velocity, and An, urn is the 
radial distance over which most of the energy is released. The ratio 
of nuclear to viscous heating can thus be written 
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where we have used equation I ll2t . neglecting terms oc (H/r) 2 . 

Due to the temperature sensitivity of nuclear burning rates, 
T = ''bum generally occurs within a factor < 2 of a fixed tempera- 
ture T = T bum , which depends primarily just on the element under 
consideration. Because the radius and temperature of burning are 
related by 



Figure 7. Bolometric light curve of supernova-like emission from WD- 
NS/BH mergers, powered by the decay of 56 Ni in the wind-driven ejecta. 
Three models shown correspond to the merger of [1] a 0.6M B C-O WD 
with a 1.4M NS (solid; v w a 2.8 x 10 4 km s" 1 ; M CJ a 0.5M o ; M Ni = 
6 X 10~ 3 Af Q ; Fig. [JJ; [2] a 0.3M o He WD with a 1.2M Q NS {dotted; 
v w a 2.2 X 10 4 km s~ ! ; M ej =s 0.3M G ; M Ni = 3 X 10~ 4 Af o ; Fig.fJ); and 



[3] a 1.2M Q O-Ne WD with a 3M Q BH (dashed; 
■■ 5 x 10- 2 M o ;Fig.[5) 
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in gas and radiation-dominated regimes, respectively, we can write 



*/visc 



H 

Art™- 



independent of M or the mass of the central object. However, this 
result holds only if gas pressure dominates (M > lCr 3 (ce/0. 1)M Q 
s" 1 ; eq. 1241 ). At lower accretion rates, radiation pressure domi- 
nates, in which case qmc/9-nsc decreases oc M° A M~ 0A . 



5 ELECTROMAGNETIC COUNTERPARTS 

In this section we use our results from ij4] to assess possible elec- 
tromagnetic (EM) counterparts of WD-NS/BH mergers. 
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where we have assumed that H/r = 0.4 and a mean molecular 
weight /i = 2. Assuming nominal gas pressure dominance, the ratio 
of radiation to gas pressure is given by 
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Turbulent mixing likely sets a lower limit of Ar bum > H on 
the radial extent of burning. On the other hand, burning is unlikely 



to occur over a much larger radial region Arb um » r bul 



2H 



due to the sensitive temperature dependence of the nuclear reaction 
rates and the power-law behavior of T(r). The bracketed quantity 
in equation J24t is thus a reasonable estimate of q mc /q VKC . The ap- 
proximate equality ^ nuc > ^ v i sc may be considered a loose defini- 
tion for whether a disk is in the NuDAF regime on a radial scale 

* ~ 'burn- 

As an example, consider carbon burning (Q/A & 0.38 MeV 
nucleon -1 ; T bum ss 1.5 x 10 9 K) in a NS disk with a mass fraction 
X\ 2 = 0.5, for which A6 nuc ss 0.2 MeV. Equation {24]! shows that 
<?nuc < qvix for M = 10 _3 M Q s - ' and a = 0.1, a result seemingly 



(24^ previously discussed EM counterpart of WD-NS/BH mergers is 
a long durat ion Gamma -Ray Burst, p owered by accre tion onto the 
central BH dFrver et al.|[l999l) or NS dKing et al.ll2007h . Although 
the formation of a relativistic jet from the inner disk is certainly 
possible, our calculations show that the accreted mass reaching 
the NS surface or BH event horizon may be a factor ~ 10 - 100 
smaller than that of the original WD (although see 94.4t , Among 
other things, this suggests that WD-NS mergers resulting in BH 
formation may be relatively rare. 

Even if sufficient mass accretes to produce a powerful jet, only 
a small fraction of GRB jets are pointed towards the Earth. For off- 
axis events, the prompt and afterglow emission are much dimmer 
due to relativistic debeaming. Thus, although it would be unsur- 
prising if WD-NS/BH mergers were accompanied by non-thermal 
jetted emission at some level, bright high energy emission may not 
be a ubiquitous feature. 



5.2 Radioactively-Powered Optical Transient 

A promising source of isotropic EM emission from WD-NS/BH 
mergers is a supernova(SN)-like transient, powered by the radioac- 
tive decay of 56 Ni from the wind-driven ejecta. We argued in 94.21 
that outflows from the disk form a singular homologous 'shell' on 
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large scales. Depending on the binary parameters and the wind ef- 
ficiency r] w , the ejecta has a characteristic mass M c j ~ 0.3 - 1M , 
mean velocity v e j w l-5x 10 4 km s -1 , and Ni mass M Ni » X 56 M c j ~ 
10~ 45 - 10~ 25 M o (Table|2](. More massive WDs generally produce 
more massive ejecta, with a higher speed and larger Ni yield. This 
conclusion should be taken with caution, however; thermal insta- 
bility may alter the disk evolution, especially for high mass (and 
possibly hybrid) WDs (j g3) , 

Figure [7] shows bolometric light curves of Ni decay-powered 
tr ansients from W D-NS/B H mergers, calcul ated using the model 
of lKulkarnil feOOSt) and lMetzger e t al. (2008a) for ejecta properties 
corresponding to different examples of WD-NS/BH systems from 
ij4] In all cases the light curve peaks on a timescale ~ week. By 
contrast, the peak luminosity L\, i varies from ~ 10 39 ergs s _I to 
~ 10 415 ergs s -1 (absolute magnitude Mb - -9 to = -15), as re- 
sults primarily from the large spread in 56 Ni mass, which depends 
primarily on the mass of the disrupted WD. 



5.2.1 Comparison to observed subluminous Type I SN 

At least in the case of C-O and O-Ne mergers, both the brightness 
and duration of the predicted transient are broadly consistent with 
the properties of recently discovered subluminous Type I SNe such 
as 2005E, 2008ha, and 2010X (JTJ. On closer inspection, however, 
connecting individual events (or classes of events) to WD-NS/BH 
mergers results in several potential difficulties. 

Although the ejecta mass, Ni mass, and presence of C, O, 
Si, S, Ca, Fe from C-O/O-Ne WD mergers are all broadly con- 
sistent with the observed pro perties of SN 2008ha ( Valen ti et al.l 
l2009l : lFolev et alj|2009ll201fj|) . the velocity of the ejecta < 10 4 km 
s -1 inferred for this and other '2002cx-like' events are much lower 
than the predictions of our baseline models. One caveat is that our 
calculations of v c j assume that winds from all radii in the disk con- 
tribute to a single homologous body of ejecta, whereas in reality v c j 
could be smaller if the fast wind from the inner disk escapes along 
the pole without efficiently coupling its kinetic energy. 

The ejecta mass, Ni mass, velocity, and presence of He, Ca, 
and O of our C-0 and He WD models are similarly consistent with 
the properties of SN 2005E. However, the large Ca mass > 0.1 Mp 
inferred in the ejecta of this and related events jPerets et al.112010]) 
is much higher than we predict. If this disagreement can be rec- 
onciled, we speculate that the mysterious location of 05E-like ob- 
jects in the outskirts of their host galaxies could be explained by 
our model if WD-NS binaries are given a 'kick' during their super- 
nova, which removes them from the disk of the galaxy by the time 
of merger. 

For SN 2010X, again many of the elements produced by 
WD-NS/BH mergers are seen in the spectra (and some, like C 
and other interme diate mass elements, are not expected in alter- 
native .la models; iKasliwal et al. 2010). However, Ti is also ob- 
served, despite its low predicted quantity in our models. Further- 
more, n one of our models pr oduce enough 56 Ni to explain SN 
2002bj JPoznanski et"aT1l2010l) . an event which otherwise shares 
several properties with 2010X. Again, these conclusions must be 
moderated due to our ignorance of the outcome of high mass ratio 
mergers and the uncertain effects of thermal instabilities. 

As a final note, although this possibility has been neglected 
thus far, free nucleons from the very inner disk may also con- 
tribute a small fraction of the ejected mass (Fig.[T}. These may be 
lost to winds heated by viscous dissipation (as discussed here) or 
by neutrino irradiation from the very inner disk or boundary layer 



(e.g. lMetzger et alj20o"8l . lMetzger et al.l2008bT) . Once neutrons de- 
cay, the net effect is contamination of the ejecta with hydrogen. 
Since even a small quantity of H may be detectable due to its strong 
lines, its presence in an event otherwise classified 'Type I' would 
support the WD-NS/BH merger model, due to its unlikely presence 
in other WD models (its absence would, however, not rule out the 
model). 



5.3 Radio Transients 

A final source of transient emission from WD-NS/BH mergers is 
non-thermal synchrotron emission powered by the deceleration of 
the ejecta with the surrounding interstellar medium (ISM). Shock 
emission peaks once the ejecta sweeps up its own mass in the ISM; 
this occurs on a timescale 

^H^kTii^f^Y 13 ^ (26) 

where n is the ISM density. Following Nakar & Pirarj d201l[) (their 
equation [14]), we estimate that the number of radio transients de- 
tectable with a single 1.4 GHz snapshot of the whole sky down to 
a limiting flux F tim is given by 

f_5_ v*. y 3 (±.r (_^v 15 f * ) (27 ) 

\cm- 3 /\0.l/ \0.l/ \0.1mJyj ^lO^yr" 1 / 

where R ~ 10~ 5 - 10~ 3 yr~' is the merger rate per galaxy (see 5j2} 
and we have assumed that relativistic electrons are accelerated into 
a power-law energy distribution with index p = 2.5. Here e e (e B ) are 
the fraction of the energy density behind the shock imparted to rel- 
ativistic electrons and magnetic fields, resp ectively, norma lized to 
values characteristic of radio supernovae (IChevalierl 1 1 9981) . Equa- 
tions J26t and Ml\ sh ow that future w ide-field radio surveys at ~ 
GHz frequencies (e.g. lBower et al.|[2010t) will preferentially detect 
those events with the highest ejecta mass and velocity, as charac- 
terize higher mass WD mergers in our model. This is due both to 
the sensitive dependence of N\a on M c j and v e j and the requirement 
that t ieQ be sufficiently short to produce an appreciable change in 
the radio brightness over the characteristic lifetime of the survey. 



6 DISCUSSION AND CONCLUSIONS 

Although the effects of nuclear burni ng on accretion have 
been explored in prev ious work (e.g. iTaam & Frvxelll 1 19851 
Ichakrabarti et al.l 119871) . these efforts focused on systems with 
much lower accretion rates, such as X-ray binaries. Unphysically 
low values of the disk viscosity a < 10~ 10 are necessary in these 
cases to achieve sufficiently high densities for appreciable nucle- 
osynthesis. Here we have shown that nuclear burning has an im- 
portant effect on the dynamics of accretion following the tidal dis- 
ruption of a WD by a NS or stellar mass BH, even for more re- 
alistic values of a > 0.01 - 0.1. Under conditions when the heat- 
ing from nuclear burning is comparable to, or greater than that, re- 
leased by viscous dissipation, we have introduced the concept of a 
'Nuclear-Dominated Accretion Flow' (NuDAF; see the discussion 
surrounding eq. 11241 for a more concrete definition of this regime). 
The dependence of the disk thermodynamics on nuclear burning is 
particularly acute because inflow is already radiatively-inefficient 
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and hence marginally bound jNaravan&Yllll994l) . even without 
an additional heat source. 

In j}5]we discussed several sources of transient emission as- 
sociated with WD-NS/BH mergers. In addition to their EM signa- 
tures, WD-NS/BH binaries are considered a promising source of 
gravitational wave emission at frequencies <mHz. Estimates sug- 
gest that a spaced-based interferometer such as LIS A could resolve 
~ 1 - 100 WD-NS/BH systems within our g alaxy jNelemans et al.l 
1200 it ICooravl 120041; IPaschalidis e7ai1l2009» . Unfortunately, most 
of the discovered systems will require > 10 4 years to merge. Al- 
though the odds are not favorable, if a WD-NS/BH binary in the 
Milky Way or a nearby galaxy were sufficiently compact to merge 
on a shorter timescale, it should be relatively easy to detect and 
localize. EM counterpart searches could then be triggered at the 
expected merger time. 

The results presented here may be applied to other con- 
texts. One immediate possibility is the tidal disruption of a 
WD by an intermediate m ass BH (e.g. iRosswog et ail 120091 ; 
IClausen & Eracleousl 1201 lh . For nearly circular orbits prior to 
merger (as we have discussed), burning is unlikely to be impor- 
tant if the primary mass M is too large. This is because when 
radiation-pressure dominates (as is likely in this case; eq. 1251 ). 
the ratio of nuclear to viscous heating g nU c/?visc <* M -0 8 (eq. 1241 ) 
is small due to the deep potential well of the BH. On the other 
hand, for highly elliptical orbits with small pericenter radii, tidal 
forces during the disruptio n can themselves trigger a nuclear run- 
away dRosswog et aTT2 009 ) ; because the initial orbit has a very low 
binding energy, the nuclear energy released will have a more sig- 
nificant effect on what matter is unbound. 

A perhaps more promising application of our results is to ac- 
cretion following the core collapse of a massive star, as in the 
'colla psar' model for long duration GRBs (MacFadven & Wooslevl 
1999). Notably, the outer layers of the Wolf-Rayet progeni- 
tors of long GRBs are pred icted to have a He-C-0 composi- 
tion jWooslev&HegeJl2006h . similar to that of a 'hybrid' WD 
(Fig. (4))- Due to its low temperature threshold, the reaction 
4 He+ 16 0— » 20 Ne+y may have a particularly important influence 
on the disk, if material circularizes at a sufficiently large radius. 
Note that current numerical simulations of collapsar disks gen- 
erally neglect the effects of nuclea r reactions, except perhaps at 
very high temperature > 4 x 10 9 K fMacFadven & Wooslevll 19991 ; 
iLindner etailboiOl ; iMilosavlievic et alj|20ld) . Our results show 
that may not be valid, as nuclear energy generation can be im- 
portant, even for temperatures as low as ~ 10 9 K. The effects of 
nuclear burning may be similarly important in ac cretion following 
the m erger of a helium star with a NS or BH ( Fryer & Wooslevl 
1998). The recent GRB 101225 A was interpreted as re sulting from 
accre tion following following NS-helium star merger (IThQne et al.l 
1201 lh : we note that the small Ni mass inferred from the weak SN 
associated with this event is consistent with our estimated yield 
from He WD-NS mergers. 

There is much room for improvement on the model pre- 
sented here in future work. For one, we have assumed a time- 
steady, height-integrated disk model. In reality, the disk may have 
a complex vertical structure, due in part to the interplay between 
MRI-driven turbulence and convection, the latter driven by the 
strong temperature dependence of nuclear heating. Complex time- 
dependent behavior is also expected, associated with both the secu- 
lar viscous evolution of the disk and, potentially, shorter timescale 
variability associated with thermal instability (possibly resulting 
in limit cycle behavior; 34.3b . We have assumed a simple out- 
flow model with two uncertain parameters, rj„ and Be d , whose val- 



ues depend on the mechanism and heating source responsible for 
driving outflows from the disk corona. The pressure dependence 
of wind cooling has important implications for the thermal sta- 
bility of the disk ( 94.31 1. Our calculations also neglect the effects 
of convection on energy and angular momentum transport in the 
disk ( 94.4| l, which, if important, could substantially increase the 
accretion timescale or even alter the radial structure of the accret- 
ing envelope. Even our simple model cannot be directly applied 
to the merger of massive WDs with NSs, since nuclear burning is 
already important during the disk's formation itself. Here, multi- 
dimensional numerical simulations that include the effects of nu- 
clear reactions will be required to fully characterize the evolution 
and fate of these systems. 

Finally, the discussion of optical and radio EM counterparts in 
95.2l and j ]5.3l was necessarily limited in scope. In future work we 
will pursue more detailed calculations of the predicted optical light 
curve and spectra, as will be necessary to more precisely compare 
the predictions of our model to observations of subluminous Type 
I SNe ( 95.2. U . This work will also address the expected detection 
rates with current and upcoming optical and radio transient surveys, 
given the current constraints on the merger rates of different WD- 
NS/BH systems. 
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